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Abstract 

We study a two-dimensional superconductor in close proximity to a two- 
dimensional metallic sheet. The electrons in the superconducting sheet are 
coupled to those in the metallic sheet by the Coulomb interaction only. We 
obtain an effective phase-only action for the superconductor by integrating 
out all the electronic degrees of freedom in the problem. The Coulomb drag 
of the normal electrons in the metallic sheet is found to make the spectrum 
of phase-fluctuations in the superconductor, dissipative at long wavelengths. 
The dissipative co-efficient rj is shown to be simply related to the normal state 
conductivities of the superconducting layer (errand the metallic sheet (erg) 
by the relation r? oc a $° E . 
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The study of phase fluctuations in disordered superconductors and Josephson junction 
arrays has received a great deal of attention [1,2,3] over the years as a mechanism for the 
superconductor-insulator transition [4,5] in systems with mesoscale disorder. The inho- 
mogenities in these systems exist on a length scale much larger than the superconducting 
coherence length. The pairing mechanism is therefore believed to be relatively unaffected 
and the superconductor-insulator transition seen at low temperatures with increasing disor- 
der is believed to be caused by enhanced fluctuations in the phase of the superconducting 
order parameter. A striking feature of this transition is the connection between an increasing 
normal state resistance and enhanced phase fluctuations in the superconductor which even- 
tually destroy global phase coherence, thus leading to an insulating state. The discovery of 
a pseudo-gap phase [6,7,8] in the underdoped cuprates has given a renewed impetus to this 
field because of the possibility [9,10] that this phase corresponds to a collection of Cooper 
pairs that have lost global phase coherence. 

Recently a superconductor-insulator transition was seen by A. J. Rimberg et. al. [11] 
in a two-dimensional array of Josephson junctions held within 100 nm of a two-dimensional 
electron gas (2DEG). A novel aspect of their work was the fact that the transition was 
achieved by tuning the conductivity of the 2DEG without altering the Josephson junction 
array directly. Motivated by this experiment we study, in this paper, the phase fluctuations 
in a 2-D superconductor held in close proximity (at a distance d) from a 2DEG. We model 
the superconductor by a two-dimensional BCS model together with a random potential and 
the 2DEG as an electron gas with a random potential. We explicitly include the full long- 
range Coulomb interaction, both inplane and interplanar. The Coulomb drag caused by 
the 2DEG is shown to make the long wavelength phase fluctuations in the superconductor 
dissipative in nature. 

The dynamics of the system of electrons is described by the action 




(la) 
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where 

L sc = 5(z-d)[Y / M^r)(-^- + h s )^+ l A{ ^ T) |2 +(A(x,r)^ T (x,r)^(x,r) + /..c.)] (16), 



L e9 = 5(z)[J2xA^r)(-^- + h e ) X a(x,T)\ (lc), 

and 

8tt 

The electrons at (x, r) with spin a are represented by the Grassman field variables ^(x, r), 
Vv( x , t) and % CT (x, t), x CT (x, r) in the superconducting layer (at z = d) and the 2DEG (at 
2 = 0) respectively. Here 

-a 2 V 2 

ks = ^m ieA o(x,d,T) + F s (x) - 4 (2a) 



and 



-a 2 V 2 

h e = -z ieA (x, 0, r) + V e (x) - e e F (26). 

2m 2 



Thus, L sc includes the electronic kinetic energy and the coupling of the superconducting 
electrons at z = d to the Coulomb potential (A ) as well as a random potential (V s ). 
The field A is the auxilliary Hubbard- Stratonovich field obtained from the BCS contact 
interaction and g is the strength of the attractive interaction. L e9 describes the 2DEG at 
z = together with its coupling to a random potential (V e ) and the Coulomb potential. 
Iff gives the electric field energy of the system. We do not consider transverse vortex-like 
fluctuations and so the vector potential can be set equal to zero. 

At low temperatures we can ignore fluctuations in the amplitude of the superconducting 
order parameter and make the replacement [12, 13] A(x, r) = A exp [i9(x, r)] where A is 
the mean field value of | A(x, r) |. Then on going over to a gauge in which the order 
parameter is real [13], (i.e. making the transformation ^(x, r) — > exp [^]^ (T (x, r)) we find 
that L sc becomes 

L sc = 5{z - d^L 1 + L 2 } (3a) 
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where 



and 



Ft j + a™ ^2 

^ = E&hr + ^ - 2eA ) + U 9 + 2 J + K(x) - e|# CT (36) 
^ ar 2 zmi 

L 2 = ^ + A (^| + /i-c) (3c) 
We now proceed to integrate out the fermions. We first consider L eg . Then we have 

DxD X e~ se9 = e' s ^f [M (4a) 



'3, 



where 

S e ° = J o drj d 6 rL ea (46) 

. We then compute S^jj-[A ] to second order in A by performing a functional Taylor expan- 
sion of Slfj. There is no contribution to S e e 9 ^ at first order in A as the contribution coming 
from the electrons is exactly cancelled by the positive ionic background whose conribution 
has not been explicitly written earlier. It contributes a term of the type ieA (x, 0, r)n eg 
where n eg is the average electronic density in the 2DEG. 
At second order we find that S e ff is given by 

S% f = y £ dr J* dr' j d 2 x J d 2 x'A (x, 0, r)A (x / , 0, r')P e ° 9 (x - x', r - r') (5). 

Here 

P e ° fl (x - x', r - t') = (T T [n es (x, r)n e ,(x', r')\) (6), 

n eg (x., r) is the electron density fluctuation operator and averages are performed with respect 
to S eg after setting A to be zero. 

We now integrate out the superconducting electrons. Once again we write 

J D^Dije~ ssc = e -W'*>] (7). 

Proceeding as before, we find at first order there is no contribution to 5*1// • At second order 
we find 

S s e c ff = ±J P drdr' J rf 2 a;rfV[P s c (x-x',r-r / )/(x,r)/(x / ,r / )+P(x-x / ,r-r , )V^(x,r)-V^(x / ,r / )] 

(8). 



Here 

/(x,r)=0(x,r)-2eA)(x,d,r) (9a), 
P° (x - x', r - t') = (T r [n sc (x, r)n sc (x', r')]> (96), 

and 

£>( x _ x ', r _ r ') = I^( x _ x / )5(r _ r / } _ _L<T T [p Xl (x, rK(x', r')]> (9c). 

Tfl\ Tfl-y 

n sc represents the electron density fluctuation operator in the superconducting state and p Xl 
is the momentum density operator along the x\ direction. All averages are performed with 
respect to L sc with A and 9 set equal to zero. 

Including the electric field energy in L e f we find that the effective action for the system 
becomes 

S eff = SZj f + S% f + S e f (10a) 
where S e e 9 ^ and S*jf are defined in Eqs. (5) and (8) respectively and 

Jo 8tt y j 

Varying the effective action S e ff with respect to A (x, z,t) we obtain the equation of 
motion to be 

V 2 A 

47T 

where 



6(z-d)X(x,T) + 6(z)Y(x,T) (11a) 



X(x, t) = -± £ dr' J dVP° (x - x', r - r')( - 2eA (x', d, r')) (116) 

and 



F(x,r) = e 2 J^dr' J dVPj(x-x',T-r')A (x' 1 ) T') (11c) 

Going over to Fourier space by using the transformations 

1 . f d^q 

#( x > r) = - ^ J J7^y exp i[q • x - i/ m r]0(q, i/ m ) (12a) 

and 

^ ^ p d? q f dk ~ 

A q (x,z,t) = -J2 J J2ny J 2^ 6Xpi ' q ' X + kz ~~ u rnr]A (q, k, is m ) (126) 



we find that Eq. (11) can be rewritten as 

A (q, k, v m ) = 2 j- , 2 [-27re exp [-ikdjX^q, z/ m )-47re 2 exp [-ikd]X 2 (q, z/ m )-47re 2 X 3 (q, v m )\ 

(J T ft 

(13a). 

Here 

^i(q, ^m) = P s ° c (q, z/ m )iz/ m 6»(q, i/ m ) (136), 
X 2 {q, v m ) = P s ° c (q, z/ m )A)(q, 2 = d, ^m) (13c), 

and 

v m ) = P° g (q, u m )A (q, z = 0, i/ m ) (13d). 

Eq. (13) is an integral equation for A (q, k, is m ). Solving for A (q, d, v m ) and A (q, 0, i/ m ) 
we find, after some algebra, the relations 

A)(q, d, O [l+^P°(q, i/ m )]+^ e -*J* ( q , i/ m )4>(q, 0, i/ m ) = Z^pLp°( q , ^ m )0(q, „ m ) 

(14a) 

and 

A (q, 0, „ m )[l+^J* (q, z/m )] + ^! e -^p0 (q? ^(q, d> ^ = Z^ e -^ p o (q> ^ ) 

(146). 

Solving Eqs. (14a) and (14b) for A (q, d, u rn ) we find that 

A (q, d, i/ m ) = — — (15a) 

1 + rl2 



where 

2q q 

and 



Ri = ^ - ^^° c (q, 0^(q, 0(1 " ^ d ) (156) 



^2 = — (q, O + P e ° 5 (q, + -^PsM, 0^(q, 0(1 " e" 2 * 1 ) (15c). 

Substituting the equation of motion (Eq. (11) ) in the expression for S e ff (Eq. (10)) we 
find that S e ff becomes 

S eff = S lff + S eff ( 16 «) 
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where 



5^ = i jf rfrrfr' / d 2 ^VP° (x - x', r - r')^/(x', r') (166) 
and 

5 2 // = 1 jf^ droV y d 2 xd 2 x'D(x - x', r - r') V0(x, r) • V'0(x', r') (16c). 

For simplicity, we first consider the limits — > oo and — > 0. In the former case we find 
that 

A ° (q '"'" m) " l + (2vreVg)P°(q,^) (17) - 
Notice that in this limit the layers become decoupled and P e ° 9 drops out of the expression 

for A in the superconducting layer. Substituting Eq. (17) in the earlier relation (Eq. (16b) 

) for Sljf we get 

ol J_V- f <?Q "It I g(q^n) | 2 ^"(q^m) ( -, R) 

e// 8/3 £ri (2tt) 2 1 + (27re 2 /g)P° (q, i/ m ) 1 J 

This in combination with S 2 jj (Eq. (16c)) is the usual action [13] for phase fluctuations in 
a two-dimensional superconductor and at long wavelengths the phase fluctuations obey a 
dispersion relation proportional to g 1 / 2 . 

Assuming that the onset of superfluid order doesn't affect the electronic compressibility 
very much, we have the relation 

1 + (27re 2 /g)P s ° c (q, v m ) « e s (q, v m ) (19) 

where e s (q, v m ) is the dielectric function of the superconductor in its normal state. This 
dielectric function, in turn, is related to the normal state conductivity (as) of the supercon- 
ductor by the relation 

47rft<7 5 (q, i/ m ) Airha s (q, v m ) 
e s (q, v m ) = l + : : « : : (20a) 

I | | Vm I 

at low frequencies. In a similar manner, at low frequencies we have the relation 

o 47raag(q,z/ m ) 
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Making these substitutions in Eq. (18) we find that S e ff reduces to 

= ^ E / 70 1 "(a, o l 2 [f| + «(4 <^ 2 ] (2D 



which is the plasma action for a two-dimensional superconductor. 

We now turn our attention to the other limit of d — > 0. In this case, we find that 

A (a d u \ - -7r(e/g)i^ m Pj ) c (q, v m ) 

Using the result of Eq. (22) for A and taking the low frequency limit as before we find that 
S e ff becomes 

1 ^ f d 2 q 2 has(J E q \ v m \ 2 

where <je is the conductivity of the 2DEG and we have replaced the wavevector and frequency 
dependent conductivities by their d.c. values in the low energy limit. In writing Eq. (23) 
we have made use of relations similar to Eqs. (19), (20a) and (20b) for the electron gas as 
well. 

It is clear from Eq. (23) that in this case the phase fluctuations are dissipative in nature 
and the corresponding viscosity co-efficient (rj) is related to the conductivities of the two 
types of electrons by the relation 

V oc ■ (24). 

We will now consider the general expression for A (q, d, v m ) (Eq. (15)). In this case we 
find that S e ff becomes 

*» = ^/w' 9(q ' Vm) 'TO* + B(q ' Vm)A (25a) 

where 

ii(q, v m ) = PsM, ZAn)e e3 (q, v m ) (256) 

and 

r 2 (q, v m ) = e es (q, v m ) + e sc (q, v m )-l + (2vre 2 /g)P°(q, O^fa' 0(1 " e~ 2<zd ) (25c). 
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Here e es (q, v m ) = 1 + ^-P e ° (q, v m ) is the dielectric function of the 2DEG. 

It is clear from Eq. (25) that dr 1 defines a crossover scale and phase fluctuations with 
q d^ 1 (e~ qd — > 0) behave like the d — > oo limit considered earlier i.e. they are propogating 
modes. In this case Z(a, u m ) = ^^< Vm \ reduces to — p sc(<i^m) — ac ti n for these modes 
becomes independent of -P e ° 9 (q, f m ) showing that the short wavelength modes are unaffected 
by the presence of the metallic layer. On further taking the low frequency limit, as before, 
Z(q, u m ) becomes as in the special case of d — > 00 considered earlier. 

The opposite limit of q <C d" 1 is like the case of d — > treated earlier. In that case 
e -2<?d ]_ an( j _2'(q ) reduces to 

%,J=^^L (26). 
e es (q, i/ m ) + e sc (q, v m ) - 1 

Now taking the low frequency limit we find that Z(q, v m ) becomes 

Z(q.O= | y^'^fa^' (27). 

Further in the low energy limit we can replace the conductivities £Js(q, v m ) and a E (q,u m ) 
by their d.c. values a s and a E - Then it is clear from Eqs. (25) and (27) that the long 
wavelength modes are dissipative in nature. Their behaviour is identical with the d — > 
limit considered earlier. Notice once again that the co-efficient of dissipation (77) is related 
simply to the normal state conductivities of the superconducting sheet (as) and the metallic 
sheet (cfe) by the simple relation r\ oc ■ Thus the amount of dissipation can be tuned 

by changing the conductivity in either sheet. 

Finally, let us summarise the main results of this paper. We have obtained the effective 
action for the phase fluctuations in a two-dimensional superconductor in close proximity to 
a metallic sheet. We find that the long wavelength (q < d~ r ) modes are dissipative. The 
coefficient of dissipation is related to the normal state conductivities of both layers. The 
short wavelength (q > d~ v ) modes are unaffected by the presence of the metallic sheet. The 
wavevector d~ Y defines a crossover scale for the behaviour of the phase fluctuations to change 
from one regime (dissipative) to the other regime (propogating). 
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